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Abstract 

Dark current and multiple electron impacts (multipact- 
ing), as for example observed in radio frequency (RF) 
structures of accelerators, are usually harmful to the equip- 
ment and the beam quality. These effects need to be sup- 
pressed to guarantee efficient and stable operation. Large 
scale simulations can be used to understand causes and de- 
velop strategies to suppress these phenomenas. 

We extend OPAL, a parallel framework for charged par- 
ticle optics in accelerator structures and beam lines, with 
the necessary physics models to efficiently and precisely 
simulate multipacting phenomenas. We added a Fowler- 
Nordheim field emission model, two secondary electron 
emission models, developed by Furman-Pivi and Vaughan 
respectively, as well as efficient 3D boundary geometry 
handling capabilities. The models and their implementa- 
tion are carefully benchmark against a non-stationary mul- 
tipacting theory for the classic parallel plate geometry. A 
dedicated, parallel plate experiment is sketched. 

INTRODUCTION 

Dark current and multipacting phenomena have been ob- 
served in various accelerator RF structures, e.g. in electron 
guns, due to field emission caused by strong accelerating 
fields [1 1, and multipacting is also appearing in high-Q RF 
cavities of cyclotrons El [3). These phenomena are usually 
harmful to the equipment and beam quality, as they will 
cause galvanic etching on the surface of the cavity and thus 
cause RF breakdown. 

Multipacting in cyclotron cavities is a very disturbing 
phenomenon. The seed electrons will impact the cavity 
surface, and produce an avalanche of new electrons. Un- 
der certain conditions (material, geometry of the RF struc- 
ture, frequency and level of the electromagnetic field), the 
secondary emission yield (SEY) coefficient will be larger 
than one and lead to exponential multiplication of elec- 
trons. This kind of discharge will limit the power level, 
until the surfaces are cleaned through a conditioning pro- 
cess. However, this process is very time-consuming J2][3l. 

Large scale dark current and multipacting simulations 
based on reliable data of surface material, full size geom- 
etry model of RF structures and parallel computing capa- 
bilities, allow more thorough analysis and a deeper under- 
standing of these phenomena, even in early design stages 
of RF structures. 
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To make OPAL [4] a feasible tool for performing large 
scale dark current and multipacting simulations, we imple- 
ment a 3D particle-boundary collision test to minimise par- 
ticle searching during the tracking process. In a subsequent 
step we add surface physics models including an analytic 
Fowler-Nordheim field emission model and two secondary 
emission models, developed by Furman-Pivi and Vaughan 
respectively. The above mentioned models and their imple- 
mentation in OPAL have been benchmarked against a non- 
stationary theory J3]. A nano-second time resolved multi- 
pacting experiment is ongoing. 

MODELS 

Geometry Handling 

Testing particle-boundary collisions is crucial to both 
dark current and multipacting simulations. Since complex 
3D geometries are hard to be accurately parameterized by 
simple functions, we use triangulated surfaces, which are 
extracted from a volume mesh generated by GMSH [6|, to 
represent the complex geometries. Subsequently we can 
make use of efficient 3D line segment-triangle intersection 
(LSTI) tests to identify particle-boundary collisions. 
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Figure 1 : Schematic view of the particle-boundary early re- 
jection strategy. The dark black line represents the bound- 
ary surface, particles are coloured dots with an attached 
momenta arrow. Gray arrows visualize inward normals of 
the boundary. 

Our LSTI test algorithm is based on [7 1, and the detailed 
introduction to our implementation can be found in our pre- 
vious paper [8]. Early rejection strategies skip particles far 
away from the boundary and with equal direction of mo- 
menta and boundary normal vectors (see Figure Figure [T| 



Surface Physics Models 

Electron field emission is a major source of both dark 
current particles and primary incident particles in sec- 
ondary emission. We employ the Fowler-Nordheim (F-N) 
formula ([TJ to predict the emitted current density l9l [T0l 

where Jfn(?, t) stands for emitted electric current density 
in position r and time t. The Greek letters tp and j3 denote 
the work function of the surface material and the local field 
enhancement factor respectively. The parameter E is the 
electric field in the normal direction of surface. The pa- 
rameters A and B are empirical constants. The functions 
v{y) and t(y) representing the image charge effects ifTUl as 
a function of the Fowler-Nordheim parameter y with the 
following definition HI 

y= = 3.795 x 10"^. (2) 

V VKE If If 

In our model, we have chosen simpler approximations 
v a (y) and t a (y) for v(y) and t(y) Q): 

v a (y) = a -by 2 
t a (y) « l. 

where a and b are free parameters to fit the value of v(y). 
These approximations are valid for a large range of y, cor- 
responding to typical applied electric field ranges in RF 
guns. 

Whenever the normal components of an electric field are 
strong enough the field emission current density will be 
limited by space charge effects IfTUl . To cover this situa- 
tion we incorporated the ID Child-Langmuir law 

-mm) - 

into our field emission model. Jsc{ r ,t) denotes space 
charge limited emission current density in position r and 
time t, Eq the permittivity in vacuum, E the normal compo- 
nent of electric field on the surface and d the distance from 
the position where E is evaluated. Currently we choose d 
to be equal to the distance traveled by emitted particles in 

one time step, i.e., d — where At is simulation 

time step. In each time step, the emitted current density 
J(r, t) for each surface triangle will be the smaller one of 

J FN (r,t) and J sc (r,t), 

J(r,t) = min{J FN (r,t), J sc (r, t)}. (4) 

We implemented two secondary emission models. The 
first one is a phenomenological model developed by M. 



A. Furman and M. Pivi [11], generating various genera- 
tions of secondary electrons: true secondary, rediffused or 
backscattered. 

The other secondary emission model is based on a 
secondary emission yield formula developed by Vaughan 

mini El: 

5(E,0) = S , forv<0 (5a) 

5(E,0) = 5 max {9) ■ {v^ v ) k , forv < 3.6 (5b) 

5(E, 6) = S max {0) • 1.125/w ' 35 , for v > 3.6 (5c) 

where 

E-E 

E m ax(Q) — Eo 

k = 0.56, forv < 1, 
k = 0.25, for 1< v < 3.6, 
5 max (9) = S max (0) ■ (1 + k e e 2 /2i:), 

E max {6) - E max (0) ■ (1 + k E 2 /2n). 

The secondary emission yield value for an impacting elec- 
tron with energy E and incident angle 9 w.r.t the surface 
normal is denoted as S(E,6). Parameter kg and k F denote 
the dependence on surface roughness. Both should be as- 
signed a default value of 1 .0, which appears appropriate for 
typical dull surfaces in a working tube environment. Lower 
values down to zero or higher values, up to about 2.0, are 
only valid for specific cases lfl2l . E max (0) is the impacting 
energy when the incident angle is zero and secondary yield 
reaches its maximum value. Eq is an adjustable parameter 
to make the first crossover energy at which the secondary 
yield equals to 1 be fitted to the experiment data 1 14 1. 

Implementation within OPAL 

The above models are implemented in the object- 
oriented parallel ESPIC code OPAL|g). 

Statistical data, such as particle position, momentum and 
particle type (primaries, field emitted electrons or true sec- 
ondaries), are stored in the H5hut [ 15 1 file format. In a post 
processing step, the data can be converted into legacy VTK 
formatted files [16] and processed by standard visualisation 
tools. 

Users can customise the dark current and multipacting 
models in the MAD like input file of OPAL. The input file 
format and examples can be found in OPAL user guide [4|. 

Motivated by the fact that the number of simulation par- 
ticles may grow exponentially, spanning a range of ten or- 
ders of magnitude, a particle re-normalisation technique is 
implemented. In each electron impact event, instead of 
emitting the real number of simulation particles predicted 
by secondary emission models, this re-normalization ap- 
proach emits only one particle and a scaled charge Q s . 



Where Q s — Qintident x S, the incident particle charge 
multiplied by the SEY value 5. This approach is an accu- 
rate representation of the secondary emission model which 
can be observed in the following parallel plate benchmark- 
ing cases. 

BENCHMARKS 

Benchmark Against the None-stationary Theory 

The theory we use to benchmark the described models 
is restricted to the case of the simple plane-parallel model 
of multipactor. We consider a spatially homogeneous and 
time harmonic RF field in between and directed perpendic- 
ular to the plates i.e., 



E(M) 



-zEq sincot 



,V . . 
-z — smut 
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(6) 



as shown in Fig. [2] Electrons are assumed to oscillate be- 
tween two parallel plates separated by a distance d. The x 
and y dimensions of the plates are assumed to be infinite. 
The impact of an electron with the plates is accompanied by 
a secondary emission yield, which depends on the energy 
and angle of the primary electron and material property of 
the multipactor. 
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Figure 2: The geometry of the analytical model. 

Assuming electrons are initially generated on the surface 
of the lower parallel plate, i.e., z = at t = to according 
to Fig. [2] then the equation of motion reads 



d?z 
dt 2 



e 

— E sin cot 
m 



sin tot 

m d 



(7) 



with Vq being the peak voltage amplitude between the par- 
allel plates. 

Integrating the equation |7]) w.r.t variable t, following 
03, i e- 



dz 

substituting the initial condition — | t=to 



vo, z\ t= t = 0, and using normalised variables: v u = 
eVo/mujd, X = iod/v u , u = i>o/v u , Lota = ipo, we ob- 
tain the scaled velocity and absolute position of electrons 
before they impact as: 



d . d 

z = — — sin cut + — (u + cos fojijjt 
A A 

d d 

+ j smtpo - - (u + cosLp Q )tp . 



(8) 



If we define cot = ip, £ = loz/v u and r 
([8]l can be rewritten as: 



i^ — i^o, equation 
£(<p, (fo,u) — (u + cos^o)''" + sin^o — sin(iy9o + t). (9) 



A non-stationary statistical theory, originated from a sta- 
tionary statistic theory ifTHl . for multipactor introduced by 
Anza et al. [5 1, gives a more realistic scenario than previous 
multipacting theories, since it considers the random nature 
of the electron emission velocity and models both double 
and single surface impacts, as sketched in Fig. [3] 




Figure 3: A full scenario of electron's trajectories between 
parallel plates, including both double surface (ds) and sin- 
gle surface (ss) impacts (5). 

The basic idea of the non-stationary theory can be sum- 
marised as follows: the initial velocity u of emitted parti- 
cles is a random variable, the solution of equation Q with 
respect to time r is the joint probability that an electron re- 
leased at phase ip s impacts at the opposite wall, separated 
by A, in a transit time t. As long as we know the probabil- 
ity density function (PDF) of the initial velocity u, which 
usually is a thermal distribution, then the joint PDF can be 
derived according to the rule of change of variable in prob- 
ability theory. 

Following the definitions in J5), each plate will be de- 
noted as D and U for the boundary condition £ = (down) 
and £ = A (up), respectively. Double and single surface 
impacts with D — U or U — D and D — D or U — U tra- 
jectories will be denoted as ds and ss respectively. Other 
definitions for the non-stationary theory are given in 0, 
and listed here in Table[T]for convenience. 

Table 1: Non-stationary theory definitions. 



Impact rate (electrons/radian) in plate 
U/D at phase ip 

Emission rate (electrons/radian) in plate 
U/D at phase ip 
Number of elec trons at time ip 
Probability density that an electron start- 
ing at plate U/D, with starting phase tp, 
experiences a double/single Surface im- 
pact in a transit phase r 
Secondary emission yield of an electron 
starting at plate U/D, with starting phase 
ip, experiences a double/single surface im- 
pact in a transit phase r 



G, 



ds/ss,U/D 



°ds/ss,U/D 



The electron emission and impact rates in each plate can 
be described by the joint PDF and the secondary emission 



yield coefficient, both of which are functions of initial ve- 
locity u, initial phase cp s and time r at which particles hits 
the plates. Details of constructing the joint PDF and inte- 
grating the emission rates and impact rates can be found in 
. The number of electrons between the parallel plates at 
phase ip then can be integrated as: 

N(<p) = ^ (Cutf) + C D {y') - W) - I D tf)) dtp'. 
Jo 

(10) 

We use the secondary emission yield curve of copper 
provided in ifTTI to benchmark the Furman-Pivi model and 
use the secondary emission yield curve of silver given in 
(5) to benchmark Vaughan's secondary emission model. 

The initial particles are equally distributed between both 
plates. The velocity of emitted particles both in the the- 
ory and in the benchmark simulation follows a Maxwellian 
distribution [5 |. We have chosen different surface material 
SEY curves, voltages, gap distances, frequencies to bench- 
mark against both Furman-Pivi's model and Vaughan's 
model. Fristly we using / = 200 MHz, V = 120 V, d = 5 
mm and the material of the multipactor is copper. The 
results with Furman and Pivi's secondary emission model 
matches the theoretical model very well, as shown in Fig. 
[4] The Vaughan's model has been benchmarked with silver 
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Figure 5: Time evolution of the electron density pre- 
dicted by the theoretical model and the OPAL simulation. 
Vaughan's secondary emission model with and without re- 
normalization is used. Model parameter are: / = 1640 
MHz, V = 120 V, d = 1 mm. 



and without surface treatment within one RF cycle. The 
electron multiplication without surface treatment is 5 or- 
ders of magnitude larger than in the case with surface treat- 
ment. 
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Figure 4: Time evolution of the electron density pre- 
dicted by the theoretical model and the OPAL simulation. 
Furman-Pivi's secondary emission model with and without 
re-normalization is used. Model parameter are: / = 200 
MHz, V = 120 V, d = 5 mm. 

SEY data. Again very good agreement of the model and 
simulation can be observed in Fig. [5j 

APPLICATIONS 

Multipacting phenomenas of the CYCIAE-100 H" AVF 
cyclotron under construction at the China Institute of 
Atomic Energy (CIAE) |fl9l are investigated. 

We have chosen two different SEY curves of copper, 
with and without surface treatment, both in agreement with 
equation |5j). 

Figure [6] shows the time evolution of the electron den- 
sity. Multipacting has been observed in both cases, with 
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Figure 6: The time evolution of particle numbers in both 
with and without surface treatment cases, the phase lag of 
both cases are degree. 



For visualising the positions where multipacting hap- 
pens, OPAL will dump the impact position and current of 
incident particles into a specified file using the efficient par- 
allel H5hut file format 03]. The simulated hot spot of the 
RF cavity of CYCIAE-100 cyclotron is shown in figure|7] 




Figure 7: Simulated hot spot of RF cavity of CYCIAE-100 
cyclotron. 

The Dark current module of OPAL has also been used 
to study the Dark current of the CTF3 gun in SwissXFEL 
project ETIl . 

CONCLUSION AND OUTLOOK 

RF components with arbitrary, closed structures, like 
electron guns and cyclotron resonators can be modelled in 
OPAL to numerically study dark current and multipact- 
ing phenomena. Detailed benchmarks against the non- 
stationary multipacting theory, enables predictive numer- 
ical studies of multipacting phenomena. Post processing 
utilities, like the hot spot visualisation and time evolution 
of particle distribution, are useful to locate the multipacting 
zones. Further multipacting simulations are planed for the 
CYCIAE-100 cyclotron cavity on different RF power lev- 
els to better understand the multipacting behaviour of the 
cavity during the RF conditioning process. The aim is to 
shorten the time consuming RF conditioning process, by 
localised surface treatment (painting) predicted by OPAL 
simulations. 

Time Resolved Multipacting Experiment 

A dedicated nano-second time resolved multipacting ex- 
periment is ongoing to benchmark the models. This dedi- 
cated experiment is using a A/4 transmission line resonator, 
which works around 73 MHz. A sinusoidal electrical field 
will be established between parallel plates, in the conduc- 
tive end of the resonator. An electron pickup is mounted 
through a hole in the centre of the ground plate to collect 
the multipacting electrons. It is expected that the collected 
multipacting current will have the same pattern as predicted 
by the none-stationary theory and the OPAL model. The 
experiment is still ongoing, detailed results will be pub- 
lished in a forthcomming paper. 
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